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EVALUATION OF SPACE SHUTTLE MAIN ENGINE 
FLUID DYNAMIC FREQUENCY RESPONSE CHARACTERISTICS 


by 

T. G- Gardner 
ABSTRACT 

In order to determine the POGO stability characteristics of the Space 
Shuttle main engine (SSME) liquid oxygen (LOX) system, an evaluation 
of the fluid dynamic frequency response functions between elements in 
the SSME LOX system was performed, both analytically and experiment- 
ally. This report acquaints the reader briefly with the POGO testing 
program and, more specifically, documents, as a user note, the POGO 
data analysis software. POGO refers to the effect the dynamic inter- 
action between devices in the LOX system has on fluid/mechanical 
vibration of the system. For the experimental data evaluation, a 
software package was written for the Hewl et t-Packard 5A51C Fourier 
analyzer. The POGO analysis software consists of five separate seg- 
ments. Each segment is stored on the 5^51C disc as an individual 
program and performs its own unique function. The POGO analysis 
software includes two separate data reduction methods, a signal 
calibration, coherence or pulser signal based frequency response 
function blanking, and automatic plotting features. The 5A51C allows 
variable parameter transfer from program to program. This feature is 
used to an advantage and requires only minimal user interface during 
the data reduction process. Experimental results are also included. 
Comparison of experimental results with the analytical predictions 
permits adjustments to the general model in order to arrive at a 
realistic simulation of the POGO characteristics. 


SPACE SHUTTLE MAIN ENGINE SYSTEM DESCRIPTION 


The Space Shuttle propulsion system consists of two solid rocket boosters 
(SRBs) and three main engines (SSMEs). The SSMEs are high-performance, 
liquid-propellant, variable thrust rocket engines, operating at high 
temperatures, high pressures, and high rotational speeds. Each SSME 
operates at a chamber pressure of approximately 3000 psla to produce a 
sea level thrust of 375,000 pounds and vacuum thrust of 470,000 pounds 
and operates over a range from 50 to 109 percent of the rated power. 

Each engine consists of ten major components: two preburners, four 

turbopumps, the hot-gas manifold, the main injector, a heat exchanger, 
and the main combustion chamber. All other components are attached 
structurally to the hot-gas manifold, and the entire arrangement is 
called the SSME powerhead. 

A key to the cost effectiveness of the Space Shuttle concept is hardware 
reusability. As a result, system reliability is of paramount importance. 
The SSMEs have been subjected to extensive hot firing and flow tests. 
Under these extreme operating conditions, system failures and malfunc- 
tions have occurred due to the seif- induced dynamic environment. These 
failures have ranged from subcritical wear of component bearings to 
catastrophic explosion and fire resulting from the intense pressure 
oscillations and dynamic stresses occurring in pumps, valves, and/or 
propellent lines. The components in the liquid oxygen (L0X) system 
include the low pressure oxidizer pump (LP0P), POGO suppression system 
(or accumulator), the high pressure oxidizer pump (HP0P) , and the main 
combustion chamber (MCC). The LOX flows from the external tank ( ET ) 
through the above components and is combined with the liquid hydrogen 
(LH.,) fuel at the inlet of the MCC by the main injector. 

The IP0P is an axial-flow pump driven by a six-stage turbine and powered 
by the LOX. During engine startup and mainstage, the LPOP maintains 
sufficient pressure in the LOX line to permit the HP0P to operate at 
high speeds without an inducer and without cavitation, even under 
worst-case conditions. 


The POGO suppressor is a gas-filled accumulator, which serves as a 
capacitance in the LOX flow circuit. The unit is designed to prevent 
low-frequency oscillations, transmitted from the Space Shuttle vehicle, 
from being transmitted into the HPOP and ultimately to prevent MCC 
pressure oscillations. To provide suppression protection through the 
startup and shutdown transients, the accumulator is filled with gaseous 
helium. During normal engine operation, the heat exchanger provides 
gaseous oxygen (G0 2 ), as the compliant medium, to the accumulator. The 
system consists of a 0.6 cubic-foot accumulator, which is attached to 
the HPOP inlet duct, an internal stand pipe, helium precharge valve 
package, gaseous oxygen supply valve package, and recirculation isola- 
tion valves. The liquid level in the accumulator is controlled by the 
stand pipe, which is orificed to regulate the G0 2 overflow over the 
engine operating power level range. Excess gaseous and liquid oxygen 
are recirculated back to the LPOP inlet through the engine oxidizer 
bleed duct. 

The HPOP consists of two single-stage centrifugal pumps on a coimon 
shaft, which is directly driven by a two-stage hot-gas turbine. The 
main pump receives oxidizer from the LPOP discharge and supplies LOX at 
increased pressure to the LPOP turbine. The HPOP turbine is powered by 
hot-gas (hydrogen-rich steam). LOX enters the HPOP main pump through 
the main pump housing and flows through an inlet with a 50-50 flow split 
into a double-entry, coimon outlet impeller. Several sets of guide 
vanes direct the flow to the impeller inlets. The impeller has four 
full and four partial blades in each half. After passing through the 
impeller, the flow is redirected into the HPOP discharge by diffuser 
vanes. The HPOP is attached to the hot -gas manifold by a flange, which 
is canted ten degrees out from the engine centerline. The oxidizer then 
flows into the MCC. 

The MCC is a cylindrical structural chamber, which contains the burning 
propellant gases and is flange-attached to the hot -gas manifold. The 
MCC consists of a coolant liner, coolant inlet ana outlet manifolds, and 
a structural jacket. The coolant liner, which provides the coolant 


flow path of the MCC, contains 390 milled axial coolant channels, which 
are ported to the coolant inlet and outlet manifolds. This network pro- 
vides an up-pass circuit for the liquid hydrogen coolant. The chamber 
jacket prcvides the structural strength for the MCC and is approximately 
20 inches long. The jacket is formed in two matching halves, which are 
welded together and to the coolant inlet and outlet manifolds. A throat 
ring is also welded to the jacket at the MCC throat to provide added 
strength. The MCC throat follows a contraction from the hot-gas mani- 
fold of 2.96:1 and is expanded to a ratio of 5:1 before attaching to the 
engine nozzle. 

In order to validate system performance and ensure equipment reliability, 
the SSME and components have been and are present ly* undergoing extensive 
development and qualification tests. Testing of the engine and components 
is conducted at several NASA and contractor locations. Full-scale 
engine test firings for development and flight acceptance are performed 
on two single-engine test stands at the National Space Technology Labora- 
tories (NSTL), Bay St. Louis, Mississippi, and one stand operated by 
Rockwell International near Santa Susana, California. In addition, main 
propulsion testing is performed at NSTL on a stand designed to accommo- 
date the Shuttle main propulsion system element--the three-engine 
cluster, external tank, and orbiter systems. 

Testing is being performed on a continuing basis. The length of a given 
test is dependent on specific test objectives and may run from several 
seconds to over 800 seconds. During each test, comprehensive measure- 
ments are acquired to monitor system performance, including vibration, 
dynamic pressure and strain at critical engine locations. Several of 
these latter measurements are utilized on-line as emergency cutoff 
indicators, and all are recorded on magnetic tape for subsequent analy- 
s i s and eval uat i on. 
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SPACE SHUTTLE HAIN ENGINE POGO TEST PROGRAM 


As part of a program to determine the POGO stability characteristics of 
the Space Shuttle Main Engine (SSME) liquid oxygen (LOX) system, POGO 
tests were conducted by Rocketdyne engineers, in California, on the A3 
SSME test stand. In order to define the system fluid dynamic characteris- 
tics, frequency response functions [H(f)] between various components are 
calculated from pressure measurements taken at various locations in the 
system. To calculate the H(f)s, a known pressure signal must be applied 
to the system and the system response to this signal measured. A 
pulser was used to generate this dynamic pressure signal, with the tests 
usually including both sine sweeps and sine dwells. The first device 
located downstream of the pulser is the low pressure oxidizer pump 
(LPOP), followed by the accumulator (or suppressor), the high pressure 
oxidizer pump (HPOP), and the main combustion chamber (MCC). With the 
accumulator active in the system, the pulser signal is suppressed such 
that it does not reach the MCC with sufficient strength to be measured 
adequately. Consequently, tests were conducted with the accumulator 
active and inactive. The general system diagram below shows the rela- 
tive location of these devices. Pressure measurements were made at the 
pulser, the LPOP inlet, the HPOP inlet, and the HPOP discharge. 



The pulser sine sweeps ranged in frequency from 2 to 40 Hz. At the 
higher frequencies in the pulser sweep, the displacement of the pulser 
piston and the corresponding amplitude of the pressure pulse generated 
were very small. The amplitudes of the pressure signals further down 
the line were essentially buried in the noise level and very difficult 
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to measure. For the sine dwells, the magnitude of the pulsed signal can 
be increased, providing greater dynamic range in the data. The problems 
associated with the dwell testing are the increased time required to 
conduct the test and the fact that the frequency response functions can 
be defined only at the frequencies where the dwells are located. 

Because of this, a routine that will blank out invalid data in frequency 
response functions was included in the POGO analysis software. This 
blanking can be based on either the pulser signal or the coherence 
function, leaving only the frequency response function data at the 
frequencies where the pulser was operated. In the initial tests, the 
dwell frequencies were chosen at approximately 5*Hz intervals. In 
subsequent tests, the dwells were chosen closer together in order to 
better define the system frequency response. 
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POGO DATA ANALYSIS SOFTWARE 


For organizational purposes and because of total size, the POGO data 
analysis software is divided into five segments. Each segment is stored 
separately on the 5^5IC disc and runs on the 5^*51 C as an individual 
program. The 5**51C Fourier operating system allows nesting of as many 
as ten programs. This is the feature that allows these five programs to 
operate together and perform the total POGO data analysis. The software 
includes two frequency response function (H(f)] calculation routines. 
These two data reduction methods are referred to as the RATIO and DIRECT 
methods for calculating frequency response functions. In the RATIO 
technique, all H(f)s are Initially referenced to the pulser. These 
H ( f ) s are then divided, with the resulting quotient being the H(f) 
across various combinations of devices. For example, the H(f) between 
devices 2 and 3 in the following schematic may be calculated as the 
RATIO of the H(f) between devices I and 3 end the H(f) between devices I 
and 2. 



H(f). 


/: 


H (f ) 

hTF7 


3/1 


where H(f) n /| ■ - and n • output, 1 • input (I, 21. Here 

cross spectrum between the input and output, end Gj| is the 
of the input (pulser). 


G„. Is the 
n i 

autospectrum 


NOTE: The POGO software outputs the auto (power) spectra and 

cross (power) spectra as functions normalized to the bandwidth (Af) 
or as power spectral densities (PSDs) and cross-power spectral 
densities (XPSDs). However, the 5ASIC calculates and uses the 
functions in their nonnormal ized form to calculate frequency 
response functions JH ( f ) s ] and coherence functions ly 2 (f)s) [2]. 

In the remainder of this document, all spectral functions will be 
referred to as PSDs and XPSOs regardless of their actual state of 
normal izat ion. The reader should understand that all H(f) and 
Y*(f) calculations are made prior to normalization. 


7 




To calculate the y 2 (f) that is compatible with the RATIO H(f), the 
following formula is used [3]: 




3/1 Y ‘^2/l 


where y 2 (f) n/r , • 

21 . 


G_. I 2 



1 1 nn 


•, and n and I are output and input, as above [1 , 


To arrive at and Y *^3/2 u#in 9 th * d,r#ct "••thod involves 

calculating the frequency response function as the XPSD, between devices 
3 and 2, and divided by the PSD of device 2. The coherence is calcu- 
lated as the magnitude squared of the XPSD divided by both the input and 
output PSDs. These equations are as follows: 


H ( f ) 


3/2 


r(0 


3/2 



5 — rr 

22 33 


where “ XPSD, G 22 * input PSD, and G^ ■ output PSD. 


The POGO software also calculates frequency response functions, both 
RATIO and DIRECT, which have been blanked, based on either the coherence 
or the pulser signal. Blanking means that the frequency response func- 
tion magnitude and phase are set to zero at those frequencies where the 
value of the coherence, or pulser signal, does not exceed the user set 
minimum. 


As previously stated, the POGO analysis software consists of five seg- 
ments, each of which is stored on the 5^51C disc in a separate keyboard 
£rogram Hie (KPF) record. They are currently stored in successive 
records (1 - 5) on the Marshall Space Flight Center, Systems Dynamics 
Lab's 5**5IC Fourier software disc. 
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KEYBOARD PROGRAM FILE, RECORD I 

KPFI contains part I of thraa parts, which rerform tha RATIO calcula- 
tions. Part I sots up tha data block sl 2 a, calibration factors, total 
amount of data to ba takan, and dalta f (if) [fraquancy resolution or 
bandwidth]. This information is usad to calculate both tha RATIO and 
DIRECT H(f)s. Tha 5*51C has a block arithmetic mode, which allows any 
block of data to ba multiplied or divided by any number. This number, 
however, must ba an integer. If this number happens to ba a floating- 
point variable garameter (VP), tha 5**5IC simply truncates any part of 
tha VP that lias to tha right of tha decimal. In reference to this 
problem, tha POCO analysis software first multiplies each calibration 
factor by 100, then calibrates the appropriate data block, and finally 
divides the data block by 100. This allows for two decimal places in 
each calibration factor. The same typo procedure is used with regard to 
the VP that holds the Af parameter, to allow for one digit to the right 
of the decimal when calculating PSDs. After the initial setup, KPFI 
reads the analog data from the tape and stores it on the ADC Throughput 
File, beginning at record 0. KPFI then retrieves this time domain data, 
two channels at a time, and calculates PSDs for each channel. KPFI also 
calculates H(f), y 2 (f), and XPSDs for the three pairs, LPOP/pulaer, 
HPOP/pulscr, and HCC/pulser. This data is stored in data blocks 12 
through 27. Program control is then transferred to KPF2, which is part 
2 of the P0G0 analysis software. 

KEYBOARD PROGRAM FILE. RECORD 2 

KPF2 begins the calculations for the RATIO H(f)s and y-(f)s. The 5^51 C 
has preprogrammed keys that will perform the arithmetic necessary to 
make these calculations. The keys, however, perform what will be 
referred to as block arithmetic. The 5^SIC stores data in blocks inter- 
nally, in a format Hewlett-Packard calls f loatlng-polnt-by-bloc... This 
means that the data block is actually a group of integers, all with a 
common scale factor. Because of this method of data handling, the 
arithmetic operation performed by the block command is not floating 
point, but is integer arithmetic. These operations involved ! .i the 


RATIO calculations are limited to a 40-dB dynamic range. When inverting 
a number, as that number approaches zero, the inverse approaches infinity. 
It becomes apparent that if there are any numbers whose value is less 
than 0.01 in a data block, any number greater than or equal to I will be 
more than kQ dB down when that block is inverted. The 5^*51 C automatic- 
ally scales up the block scale factor to eliminate overflowing, which 
causes any number more than 40 dB down to be set to zero. Thts is a 
serious problem when doing the type of calculation involved in arriving 
at the RATIO H(f)s and y 2 (f)s. In both of these calculations, espe- 
cially the y 2 (f) calculation, the quotient in the formulations is domi- 
nated by the small numbers in the denominator. These small numbers are 
not reliable data points, and the valid data is lost below the 40-dB 
dynamic range. 

The problem was solved in this program by using a floating-point arith- 
metic scheme in the y 2 (f) calculations and a combination of floating- 
point arithmetic and a filtering technique in the H(f) calculations. 

The floating-point scheme involves extracting each channel (or data 
point) point-by-point from a data block, putting the value into a 
floating-point format (via a floating-point variable parameter), exe- 
cuting the arithmetic operation, and finally putting the datum back into 
the data block. A magnitude filtering technique is used in the H(f) 
calculation to eliminate a loss of data after the values are put back 
into the data block. The real and imaginary parts of the H(f) are 
filtered so that any values outside the range of +7 are set to zero. 

This allows a polar magnitude maximum of 10. The 40-dB dynamic range 
then will be 10~* to 10*. This filtering technique is not necessary in 
the y 2 (f) calculations because of the final quantity inversion in the 
formulation. All the values to be inverted, however, are scanned for 
values of identically zero before inversion. This is to eliminate any 
numbers actually going to the upper limit (infinity). 

The program loops through the above operation enough times to calculate 
the RATIO H(f) and y 2 (f) for each channel in the data block. This loop 
is nested inside a loop, which allows for all three sets of ratio 


10 


calculations to be computed. Part 2 of the POGO analysis software then 
transfers control to part 3. 

KEYBOARD PROGRAM FILE. RECORD 3 

KPF3 defines the variable parameter, corresponding to the ratio calcula- 
tions, necessary to operate the plot routine. KPF3 also performs the 
calculations Involved in the H(f) blanking routine. KPF3 has several 
entry points, which allow various operations to be performed outside the 
normal program flow. KPF2 flows into the beginning of KPF3 t which jumps 
to the blanking routine. The blanking routine begins by prompting the 
user to establish whether the blanking will be based upon the y 2 (f) or 
the pulser PSD. The user will enter either a 2ero, to establish coher- 
ence blanking, or a two, to define pulser blanking. This number (0 or 
2) will also be used to establish which labels will be given to the 
plots of the blanked H(f)s during the plot routine. If pulser blanking 
is chosen, KPF3 will display the pulser PSD to allow the user an oppor- 
tunity to view the signal and decide on the minimum value used for the 
blanking routine. If y 2 (f) blanking is chosen, KPF3 skips directly to 
the next step, which is a prompt, asking for the minimum value in the 
blanking routine. Next, KPF3 sets up the block that will be used for 
generating the blanked H(f). KPF3 then enters the blanking loop. This 
loop operates individually on each channel in the frequency domain data 
block. In this loop, KPF3 goes to either the pulser or y 2 (f) block and 
gets one channel at a time, beginning with channel one (channel zero is 
the first channel) and ending with the next to last channel in the data 
block. The value of this channel is compared with the values of the 
channel immediately preceding and the channel inmediately following. If 
the value of the channel in question is greatest, it is determined to be 
a local maximum. The value is then compared with the minimum value 
specified earlier by the user. If the value of this channel is greater 
in all three tests, it is chosen as a channel of interest. The corre- 
sponding channel in the H(f) data block is gathered into a complex vari- 
able parameter and stored in the data block that was prepared previously 
for the blanked H(f). When the loop is complete, the blanked H(f) will 
consist of calculated H(f) data only at those frequencies where the 
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pulser or y 2 (f) passed all the necessary tests. The remainder of the 
blanked H(f) block will be filled with zeros. Because the HP 5^*5 1 C 
stores frequency domain data with both real and imaginary information in 
a combined complex channel, the blanked H(f) contains blanked phase 
in f ormation as well as blankod magnitude. KPF3 then loops three times 
to generate a blanked version of each H(f), HPOP/LPOP, MCC/HPOP, and 
MCC/LPOP. KPF3 then jumps to the section that defines the variable 
parameters necessary to run the plot program. Finally, KPF3 transfers 
control to KPF5, the POGO plot routine. 

KEYBOARD PROGRAM FILE, RECORD 

KPFk is a single-part program, which calculates the H(f)s, HPOP/LPOP, 
MCC/HPOP, and MCC/LPOP by matching the outputs and inputs "directly" 
rather than ratioing H(f)s, which have all been referenced to the 
pulser. This program (KPFM uses several variable parameters, which 
must have been previously defined by running the first three files 
containing the RATIO software. These variable parameters are VP 1 ■ 
block size; VP2 - number of ADC Throughput records; VP2000 through 
VP2003 ■ calibration factors for channels A through D, respectively; 
VP200 i * * frequency resolution (±f) ; and VP2030 ■ minimum value for the 
blanking routine. KPF** also retrieves time domain data from the ADC 
Throughput File. These data must have been written previously on the 
disc (starting with record zero) via KPFI (RATIO software, part 1). The 
data may be written to the disc and the necessary variable parameters 
defined independently of the RATIO software. It is recommended, how- 
ever, that all programs be operated in sequence. The direct POGO 
analysis software ( KPF4 ) writes over the data blocks generated via the 
RATIO software. This program therefore should not be initiated until 
all the required ratio data blocks have been reproduced. These data 
blocks include all PSDs and XPSDs (ref. pulser), which are not recal- 
culated by the DIRECT software. Immediately after clearing the data 
blocks, KPF*» begins a loop that operates three times. The loop reads 
the correct data channels from the disc, Fourier transforms the data, 
calibrates the data, and calculates averaged PSDs via the special sub- 
routine of the 5^51C. KPF^ then calculates the H(f) and y-(f) from the 
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averaged PSD data. The H(f), y 2 (f), and XPSDs are then stored In the 
appropriate data blocks. (See "DIRECT POGO Program Menu," appendix C.) 

KPF4 now transfers control to KPF3 at the point where the H(f) blanking 
loop begins. Before jumping to KPF3, KPF4 defines the variable param- 
eters that describe the location of the DIRECT H(F)s. Once in KPF3 the 
blanked H(f)s are calculated using the same minimum value specified by 
the user during the RATIO calculations. If the user has entered the 
POGO software In the DIRECT calculations and wishes to change the mini- 
mum blanking value, the user simply jumps to label 43 of KPF4. 

Upon returning from the blanking routine, KPF4 defines the variable 
parameters necessary to operate the POGO plot program (KPF5). KPF4 then 
transfers control to the plot routine. 

KEYBOARD PROGRAM FILE, RECORD 5 

KPF5 is entitled "POGO Plot Software." This software operates with both 
the RATIO and the DIRECT software. KPF5 uses several variable param- 
eters that must be defined by the respective software whose data is to 
be plotted. The program is organized into several separate routines, 
including a plot routine, a hard copy routine, a delay routine, and a 
routine that sets up and plots all the H(f)s in three different formats. 
The program sets up and automatically plots all the data generated by 
the respective RATIO and DIRECT programs. The KPF5 prints a data block 
menu and allows the user to replot any of the data blocks. The plot 
program is geared to run with a Tektronics 4052 graphics terminal and a 
Tektronics hard copy unit. Because of various incompatibilities between 
the 545IC and the Tektronics terminal with respect to the time required 
to transmit plot commands and the time required to execute them, a delay 
routine is incorporated in the plot program. This subroutine simply 
counts from one to 50 and returns. The routine is called in several 
locations where a time delay between plot commands is necessary. 

The plot program uses several variable parameters that establish which 
program (RATIO or DIRECT) has generated the data being plotted. The 
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plot program also labels each plot with test ID Information in the upper 
right-hand corner of the plot. This information Is stored as message 
number one in each of the text buffers one and two. Text buffer #1 is 
used for the RATIO data, while text buffer # 2 is used with the DIRECT 
data. Message 11 should be changed for each new set of test data. The 
message may be changed by executing the following set of operations from 
the 5^4 5 1 C terminal : 


Step it Command 54 5 1 C Status 

1 Y 5403 (nl) busy 

This command calls text buffer nl (1 for RATIO, 2 for DIRECT) 
into the computer core and activates the editing commands. 

2 /ROI busy 

This command tells the 54 5 1 C to replace message #1. 

3 01 busy 

(Test ID) busy 

5 /■' busy 

This sequence of commands results in the test ID being recognized 
as message 01. The /* defines the end of the message. 

6 / busy 

This command terminates the text buffer edit mode. 


To allow the user to enter the program flow for special purposes, the 
P0G0 software has many labels throughout the five KPFs. Listed below 
are the various jump commands and their uses. 

J 0 I 

This command starts the P0G0 software, initially setting up the 
calibration factors, etc, and reading the analog data. 

J 10 1 

This command enters P0G0 Software Part I after all the initial 
values have been entered and the analog data has been stored on 
disc in the ADC Throughput Pile. If this command is to be executed, 
the user must first be sure that the variable parameters that con- 
tain the values for the blocksize, calibration factors, number of 
records of data, and frequency resolution Uf) are defined cor- 
rectly and that the correct time domain data is stored on the ADC 
Throughput Pile. 


1*4 


J 20 2 


This command begins the RATIO calculations for H(f)s and y 2 (f)s. 
This command may be executed if the proper data blocks reside in 
core in the correct block numbers. The data blocks necessary are 


H (f ) 

LPOP/Pulser 

Block 12 

Y 2 (f) 

II 

Block 13 

H (f ) 

HPOP/Pu 1 ser 

Block 1 A 

Y 2 ( f ) 

II 

Block 15 

H (f ) 

MCC/Pul ser 

Block 16 

Y : (f) 

II 

Block 17 


H 30 3 

This command may be used when the user wishes to recalculate the 
blanked RATIO frequency response functions. The user will have the 
opportunity to re-establish the type of blanking [pulser or y 2 (f)] 
and/or the minimum value for that blanking. The RATIO autoplot 
parameters are defined and the program control is transferred to 
the plot routine. 

J 31 3 

This command is preliminary to the P0G0 Plot Routine. This command 
enters KPF3 and defines the variable parameters necessary to run 
the automatic plot portion of the plot routine. After these param- 
eters are defined, program control is automatically transferred 
to the plot routine. This command should be used only when the 
user wishes tc enter the autoplot routine with RATIO data blocks. 
This command might be useful if, for example, the user wants to 
change the test ID message in the RATIO text buffer after the 
P0G0 software has completed the data calculations. The user can 
stop program execution and edit the text buffer, then return to 
the autoplot sequence. 

J 40 A 

This command begins the DIRECT P0G0 Software. This should be used 
only after the RATIO program has been run, defining the appro- 
priate variable parameters, and writing the time domain data on 
the ADC Throughput File. 

J A 3 A 

This command enters KPFA where necessary variable parameters are 
defined just prior to transferring control to the H(f) blanking 
routine in KPF3. This command will allow the user to calculate 
DIRECT blanked H(f)s. 
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J 44 4 


This command causes the computer to begin execution in KPF4 at 
the point where the variable parameters that control the auto- 
plot routine for DIRECT data blocks are defined. Program control 
is automatically transferred to the plot routine. 

J 53 5 

This command will cause the plot program to display on the termi- 
nal the appropriate menu corresponding to the data stored in the 
data blocks. The plot program prompts the user to set the display 
scale as desired, then press the continue button to allow the pro- 
gram to plot the data block on the graphics terminal. This com- 
mand is useful when the user wishes to return to the plot program 
after the machine has been idled and the plot routine pointers 
have been changed. 


Appendix A describes the flow of the POGO software. Appendix B contains 
a complete list of all the program steps along with a description of 
each. Appendix C contains a listing of the text record messages used 
with the POGO software. Appendix D contains a listing of both text 
buffers used with the POGO software. Appendix E contains sample data 
output from the POGO software. 
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APPENDIX A 

POGO SOFTWARE FLOWCHART 

RATIO POGO SOFTWARE - PART 1 FLOWCHART 

Step 1 \ 

vStart - Label 0 ) 


Steps 5-36 
Program Setup 


Steps 40 - 53 

Display the prompt and allow 
the user to establish the 
block size (VP1) 



Steps 57 - 61 

Di splay 
user to 
factors 

the prompt and allow the 
input the cal ibration 
(VP2000 - VP2003) 


Steps 67 - 88 

Multiply each cal factor by 100 to 
allow for two significant digits 


Steps 95 - 99 

Display the prompt and allow the 
user to input the number of data 
records in the complete session 
(VP2) 


Steps 104 - 113 

Display the prompt and allow the user 
to input Af, and multiply it by 10 
for one significant decimal (VP2004) 




RATIO POGO SOFTWARE - PART 1 FLOWCHART (Continued) 


Steps 120 - 125 

Display the prompt - "Start the 
data tape and press CONTINUE" 


''step 129 s 
\ Pauses 


Steps 132 - 136 

Initialize the ADC Throughput and 
acquire VP2 records of four channels 
each 


Step 142 
Label 10 


Steps 146 - 176 

Set up the variable parameters 
for the data handling loops 


Step 182 
Label 11 


Steps 186 - 213 
Set up the H(f) loop 


Step 219 

Calculate the number of the VP which 
contains tne appropriate cal factor 


8 









RATIO POGO SOFTWARE - PART 1 FLOWCHART (Continued) 


Step 226 

Initialize the off-line BSFA 


Step 233 
Label 12 


Steps 237 - 266 

Get the analog data from the disc, 
convert It to frequency domain, 
calibrate it, and leave the input 
in Block 1 and output in Block 2. 


Step 271 

Calculate averaged power spectrum 


Steps 283 - 307 

Calculate H(f) and -,-(f) and store 
them as follows: H(f) in Blocks 12, 

14, 16; >*(f) in Blocks 13, 15, 17 


Steps 311 - 339 

Loop to calculate PSDs and store 
them in Blocks 18 through 27 


Steps 345 - 371 
Incremental counters 
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POGO SOFTWARE BLANKING ROUTINE 

✓'"'Step 

( Label 30 } 


Steps 5-23 

Set up variable parameters for ratio 
H(t') blanking; jump to label 32 


Label 31 


Steps 31 - 91 


Set up variable parameters to plot 
ratio data with POGO plot routine 


Step 97 

Jump to label 50 in KPF 5 (plot routine) 


Start Blanking 
Routine 


Steps 106 - 115 

Display the prompt and allow the user 
to determine the type of blanking 
routine (coherence or pulser based) 


m 







POGO SOFTWARE BLANKING ROUTINE (Continued) 



Step 128 

Define the VP31 for pulser blanking 



f 

Steps 1 

Display the correct 
the user to set the 
val ue 

35 - 143 

data block and ask 
minimum blanking 


Steps 148 - 155 

Define VP39 for blanking routine 



Step 164 

Convert the block to be blanked to 
rectangular coordinates 


Steps 1 68 - 179 

Clear a block and establish qualifiers 
for the blanked H( f 1 


Default VP31 to 
coherence blanking 


II 


26 




















DIRECT PQGO SOFTWARE FLOWCHART 

/'"’’’steps 1 ■ 9\ 

V Start Label 40 J 


Steps 13-55 

Set up variable parameters 
for program operation 


Sta 

^(f) calculation loop 
Label 41 ^ 


Steps 65 - 88 

Clear the data blocks prior 
to the calculations 


Step 92 

Initialize the off-line BSFA software 


Steps 99 - 106 

Calculate the number of the 
variable parameter contain- 
ing the correct calibration 
factor 


"step 113 

■'"N 

Start - 


Data acquisition 

loop 

s "‘*>*^J.abel 42_ 

x 





DIRECT POGO SOFTWARE FLOWCHART Continued) 


Steps 117 - 146 

Read the analog data from the disc. 
Convert the data to frequency domain, 
calibrate it, and store the Input in 
Block 1 and the output in Block 2 


Step 151 

Calculate averaged power spectrums 


Step 157 

Loop to get the complete session of data 


Step 163 

Calculate H(f) and y 2 ( f ) 


Steps 169 - 187 

Store the H(f)s in Blocks 8, 10, and 
12, and the y 2 ( f ) s in Blocks 9, 11, 
and 13. 


Steps 191 - 229 

Calculate XPSDs and store them in 
Blocks 14 through 19. 


Steps 233 - 239 
Increment the block counters 








A1 5 


0 - 


0 . 


DIRECT P060 SOFTWARE FLOWCHART (Concluded) 


f 


Step 245 

Loop for all three sets of H(f) 
calculations 






;tep 251 
Label 43 



Steps 282 - 342 

Set up the variable parameters for the 
POGO plot routine 


-0 




Steps 255 - 267 

Set up variable parameters for DIRECT 
blanking routine 




Step 273 

Jump to label 32 KPF3 (blanking routine) 






0 


Step 348 

Jump to label 50 in keyboard file 5 



31 



POGO SOFTWARE PLOT ROUTINE 



1 N 
Start - Label 50 



Step 20 


Start - Plot Loop 
v Label 51 > 
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APPENDIX B 

POGO SOFTWARE LISTING 
KEYBOARD FILE 1 - RATIO POGO SOFTWARE PART 1 


1 L 

0 



5 BS 

4096 



9 CL 




12 CL 

1 



16 CL 

2 



20 CL 

3 



24 CL 

4 



28 CL 

5 



32 CL 

6 



36 Y 

5814 



40 MS 

34 



44 MS 

14 



48 Y R 

1 



53 BS 

ID 



57 MS 

14 



61 Y R 

2000 

2003 


67 Y * 

2000 

2000D 

100 

74 Y * 

2001 

2001 D 

100 

81 Y * 

2002 

2002D 

100 

88 Y * 

2003 

200 3D 

100 

95 MS 

14 



99 Y R 

2 



104 MS 

14 



108 Y R 

2004 




Label 0 

Set blockslze to 4096 


Clear blocks 0 through 6 


Call user program 5414 (blank 
terminal screen) 

Set text record pointer to zero 
(default) 

Read next text record (message #1) 

Input from terminal VP It 
(blocksize) 

Set blocksize to VP1 

Read next text record (message #2) 

Input from terminal VP2000, 

VP2001 , VP2002, and VP2003 
(calibration factors) 

VP2000 = VP2000 * 100 

VP2001 = VP2001 * 100 

VP2002 = VP2002 * 100 

VP2003 = VP2003 * 100 

Read next text record (message #3) 

Input from terminal VP2 (number 
of data acquisition loops) 

Read next text record (message #4' 

Input from terminal VP2004 (Af) 


-This means the variable of parameter 1. 
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B2 

KEYBOARD FILE 1 - RATIO POGO SOFTWARE PART 1 (Continued) 


113 

Y * 

2004 

2004D 

10 

VP2004 » VP2004 * 10 

120 

MS 

34 

20 


Set text record pointer to 20 

125 

MS 

14 



Read next text record (message #7) 

129 

D 




Display block 0 (used as program 
pause) 

132 

MS 

32 



Set ADC Throughput File pointer 
to 0 

136 

MS 

22 

4 

2D 

Write analog data to the ADC 
throughput file. No. of 
records * 2D (4 channels each) 

142 

L 

10 



Label 10 

146 

Y _ 

3 

18 


Set VP3 to 18 

152 

Y _ 

4 

19 


Set VP4 to 19 

158 

Y _ 

5 

22 


Set VP4 to 22 

164 

Y _ 

6 

23 


Set VP6 to 23 

170 

Y _ 

7 

2 


Set VP7 to 2 

176 

Y _ 

8 

12 


Set VP8 to 12 

182 

L 

11 



Label 11 

186 

CL 





189 

CL 

1 




193 

CL 

2 




197 

CL 

3 



Clear blocks 0 through 6 

201 

CL 

4 




205 

CL 

5 




209 

CL 

6 




213 

Y _ 

9 

3 


Set VP 9 to 3 

219 

Y A+ 

10 

7D 

1999 

VP10 = VP7 + 1999 

226 

Y 

41 

0 

0 0 

Call user program 41; nl=0= 
center frequency of the band of 
interest, n2=0*width of the band 
of interest, n3 B 0=ADC throughput 
start record. (Initializes off- 
line BSFA) 

233 

L 

12 



Label 12 

237 

Y 

45 

2 

0 7D 

Call user program 45; nl=2*block 


number where the data will be 
stored, n2=0*number of Hannings 
applied to the data, n3*VP7= 
channel that will be read from 
the disc. (Off-line BSFA) 
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KEYBOARD FILE 1 - RATIO POGO SOFTWARE PART 1 (Continued) 


244 

Y 

45 

1 

0 1 

Call user program 45; nl*l, 
n2*0, n3*l 

251 

* 

2 

-100 


Multiply block 2 by VP(VP10) + 

256 

; 

2 

100 


Divide block 2 by 100 

261 

★ 

1 

20000 


Multiply block 1 by VP2000 

266 

• 

* 

1 

100 


Divide block 1 by 100 

271 

SP 

1 

2 

2 

Compute the average auto- and 
cross-power spectrums, where 
block nl'input data block and 
block nl+1 ^output, n2*2 speci- 
fies dual channel, n3*2 speci- 
fies double precision 

277 

4 

12 

2D 

0 

Loop through label 12 VP2 times 

283 

CH 

1 

2 

2 

Calculate transfer function data 
where nl 'input data block, n2*2 
specifies dual channel, n3'2 
specifies double precision 

289 

X< 

1 



Load block 1 into block 0 

293 

X> 

8D 



Store block 0 into block VP8 
[H(f )] 

297 

Y A+ 

8 



VP8 * VP8 + 1 

303 

X< 

2 



Load block 2 into block 0 

307 

x> 

8D 



Store block 0 into VP8 [y 2 (f)] 

311 

L 

13 



Label 13 

315 

: 

9D 

20040 


Divide block VP9 by VP2004 

320 

★ 

9D 

20 


Multiply block VP9 by 20 

325 

X< 

90 



Load block VP9 into block 0 

329 

X> 

-90 



Store block 0 into block VP(VP9) 

333 

Y A+ 

9 



VP9 = VP9 + 1 

339 

4 

13 

4 

0 

Loop through label 13 four times 

345 

Y A+ 

4 



VP4 » VP4 + 1 

351 

Y A+ 

5 

5D 

2 

VP5 ' VP5 + 2 

358 

Y A+ 

6 

6D 

2 

VP6 = VP6 + 2 

365 

Y A+ 

7 



VP7 = VP7 + 1 


+ The minus sign in front of the VP number means to use the value in the VP 
that is specified by the value of the VP given (i.e., the value of the VP 
specified by VP10). 
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KEYBOARD FILE 1 - RATIO P060 SOFTWARE PART 1 (Concluded) 


Bit 


371 Y A+ 8 

377 # 11 3 0 

383 J 20 2 

388 . 


VP8 ■ VP8 + 1 

Loop through label 11 three times 

Jump to label 20 In keyboard 
stack 2 

End 


kk 


KEYBOARD FILE 2 - RATIO P060 SOFTWARE PART 2 


1 

L 

20 



5 

TR 

12 



9 

TR 

14 



13 

TR 

16 



17 

Y _ 

3 

3 


23 

Y _ 

4 

6 


29 

Y _ 

5 

7 


35 

Y _ 

6 

14 


41 

Y _ 

7 

12 


47 

Y _ 

8 

16 


53 

Y _ 

9 

14 


59 

Y _ 

10 

16 


65 

Y _ 

11 

12 


71 

Y _ 

12 

13 


77 

Y _ 

13 

15 


83 

Y : 

14 

ID 

2 

90 

Y A+ 

14 



95 

Y _ 

15 

17 


101 

Y _ 

17 

13 


107 

L 

21 



111 

Y TR 

21 

-5D 

0 

118 

Y TR 

21 

3D 

1 

125 

Y TR 

21 

4D 

1 

132 

Y TR 

21 

-12D 

0 

139 

Y TR 

21 

5D 

1 

146 

Y _ 

0 

0 


152 

L 

22 



156 

Y X< 

3000 

-4D 

OD 


Label 20 

Convert blocks 12 

Convert blocks 14 

Convert blocks 16 

Set VP 3 to 3 

Set VP4 to 6 

Set VP5 to 7 

Set VP6 to 14 

Set VP7 to 12 

Set VP8 to 16 

Set VP9 to 14 

Set VP10 to 16 

Set VP11 to 12 

Set VP12 to 13 

Set VP13 to 15 

VP14 = VP1/2 

VP14 = VP14 + 1 

Set VP1 5 to 17 

Set VP17 to 13 

Label 21 

Get the block qualifiers from 
block VP(VP5) and put them in 
VP21 - VP25 

Put the block qualifiers stored 
in VP21 - VP25 into block VP3 

Put the block qualifiers into 
block VP4 

Get the block qualifiers from 
block VP ( VP1 2 ) and put them into 
VP21 - VP25 

Put the block qualifiers into 
block VP5 

Set VPO to 0 

Label 22 

Get the complex data from block 
VP ( VP4 ) , channel VPO, and put 
the values into complex VP 3000 


to rectangular 
coordinates 


45 
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KEYBOARD FILE 2 - RATIO POGO SOFTWARE PART 2 (Continued) 


163 Y X< 

3001 

-5D 

OD 

170 Y X< 

2005 

12D 

OD 

177 Y X< 

2006 

-1 2D 

OD 

184 Y IF 

2005 

0 

7 

192 Y : 

2005 

1 

2005D 

199 Y If 

2006 

0 

5 

207 Y : 

2006 

1 

2006D 

214 Y A+ 

2007 

2005D 

2006D 

221 Y A- 

2007 



226 Y : 

2007 

1 

2007D 

233 Y IF 

2007 

0 

1 

241 Y _ 

2007 

0 


247 Y X> 

2007 

5D 

OD 

254 Y * 

2009 

2900D 

2902D 

261 Y * 

2010 

2901 D 

2903D 

268 Y * 

2011 

2901D 

2902D 

275 Y * 

2012 

2900D 

2903D 

282 Y * 

2013 

2902D 

2902D 

289 Y * 

2014 

2903D 

2903D 

296 Y A+ 

2015 

201 3D 

2014D 

303 Y A+ 

2016 

2009D 

201 OD 

210 v A- 

2017 

20110 

201 2D 

317 Y : 

2904 

201 6D 

201 50 

324 Y : 

2905 

201 7D 

201 5D 

331 Y I F 

2904 

7 

1 

339 Y IF 

2904 

-7 

1 

347 Y 

2904 

0 



Get the complex data from block 
VP(VP5), channel VPO and put 
them Into complex VP3001 

Get the real data from block 
VP12, channel VPO and put It 
Into VP2005 

Get the real data from block 
VP ( VP1 2 ) , channel VPO and put 
it Into VP2006 

If VP2005 * 0, skip seven steps 

VP2005 « 1/VP2005 

If VP2006 ■ 0, skip five steps 

VP2006 « 1/VP2006 

VP2007 * VP2005 + VP2006 

VP2007 * VP2007 - 1 

VP 2007 - 1/VP2007 

If VP2007 > 0, skip one step 

Set VP2007 to 0 

Put VP2007 into the real part 
of block VP5, channel VPO 

VP2009 * VP2900 * VP2902 (ac) 

VP2010 * VP2901 * VP2903 (bd) 

VP2011 • VP2901 * VP2902 (be) 

VP2012 * VP2900 * VP2903 (ad) 

VP2013 * ( VP2902 ) 2 (c 2 ) 

VP2014 = (VP2903) 2 (d 2 ) 

VP2015 = VP2013 + VP2014 (c 2 +d 2 ) 

VP2016 « VP2009 + VP2010 (ac+bd) 

VP2017 * VP2009 - VP2012 (bc-ad) 

VP2904 = VP2016/VP2015 

l(c 2 +d 2 )J 

VP2905 = VP2017/VP2015 i p bc '^l 

l_(c 2 +d : ) 

If VP2094 > 7, skip one step 
If VP2904 >_ -7, skip one step 
Set VP2904 to 0 
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KEYBOARD FILE 2 


353 Y IF 

2905 

7 

361 Y IF 

2905 

-7 

369 Y _ 

2905 

0 

375 Y X> 

3002 

40 

382 Y A* 

0 


387 f 

22 

140 

393 Y A+ 

3 


398 Y A+ 

4 

40 

405 Y A+ 

5 

50 

412 Y A+ 

12 

12D 

419 f 

21 

3 

425 J 

30 

3 

430 . 




RATIO POGO SOFTWARE PART 2 (Concluded) 


1 2 
1 1 


OD 


0 


2 

2 

2 


If VP2905 > 7, skip one step 

If VP2905 > -7, skip one step 

Set VP2905 to 0 

Put complex VP3002 Into block 
VP4, channel VPO 

VPO ■ VPO + 1 

Loop through label 22 VP14 times 

VP3 * VP3 + 1 

VP4 * VP4 + 2 

VP 5 ■ VP5 + 2 

VPl 2 - VP 12 + 2 

Loop through label 21 three times 
Jump to label 30, keyboard stack 3 
End 
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KEYBOARD FILE 3 - RATIO POGO SOFTWARE PART 3 


1 L 

30 


5 Y _ 

31 

7 

11 Y _ 

35 

6 

17 Y _ 

36 

3 

23 J 

32 


27 L 

31 


31 Y _ 

3 

1 

37 Y _ 

4 

27 

43 Y _ 

5 

16 

49 Y _ 

6 

14 

55 Y _ 

7 

12 

61 Y _ 

8 

10 

67 Y _ 

C) 

8 

73 Y _ 

10 

6 

79 Y _ 

11 

25 

85 Y _ 

12 

24 

91 Y _ 

13 

125 

97 J 

50 

5 

102 L 

32 


106 MS 

34 

82 

111 MS 

14 


115 Y R 

30 


120 Y IF 

30 

0 

128 Y A+ 

31 

31 D 

135 X< 

31 D 


1 39 MS 

14 


143 Y R 

2030 


148 Y : 

39 

100 

155 Y A- 

39 


160 L 

33 


164 TR 

350 



1 

11 


Label 30 

Set VP31 to 7 

Set VP35 to 6 

Set VP36 to 3 

Jump to label 32 

Label 31 

Set VP3 to 1 

Set VP4 to 27 

Set VP5 to 16 

Set VP6 to 14 

Set VP7 to 12 

Set VP8 to 10 

Set VP9 to 8 

Set VP10 to 6 

Set VP11 to 25 

Set VP12 to 24 

Set VP13 to 125 

Jump to label 50 In Keyboard 
File 5 

Label 32 

Set the text record pointer to 

Read next text record (type of 
blanking routine 

Input from terminal VP30 

IF VP30>0, skip 1 step 

VP31 - VP31 + 11 

Load blanking block (VP31) 

Read the next text record 
(minimum blanking value) 

Input from terminal VP2C30 

VP39 » VP1/2 

VP39 ■ VP39 - 1 

Label 33 


2 


RATIO P060 SOFTWARE FART 3 (Concluded) 


B9 



KEYBOARD FILE 3 

- RATIO P0( 

168 CL 

360 



172 Y TR 

21 

350 

0 

179 Y TR 

21 

360 

1 

186 Y _ 

32 

0 


192 Y _ 

33 

1 


198 Y _ 

3k 

2 


20k L 

3k 



208 Y _ 

3000 

0 


21k Y X< 

2032 

310 

320 

221 Y X< 

2033 

3)0 

330 

228 Y X< 

203k 

310 

3k0 

235 Y IF 

2033 

20300 

3 -2 

2k3 Y IF 

2033 

20320 

2 -2 

251 Y IF 

2033 

203kD 

1 -2 

259 Y X< 

3000 

350 

330 

266 Y X> 

3000 

360 

330 

273 Y A+ 

32 



278 Y A* 

33 



283 Y Ao 

3k 



288 # 

3k 

390 


29k Y A+ 

35 

350 

2 

301 Y A* 

36 



306 Y AO 

31 

310 

300 

313 f 

33 

3 



319 < 
322. 


Clear block BP36 

Cot the block qualifiers from 
block VP35 and put thorn into 
VP2I - VP25 

Put tho block quo) I flora « to rod 
in VP2I-VP25 into block VP36 

Sot VP32 to 0 

Sot VP33 to I 

Sot VP3k to 2 

Label 3k 

Sot VP3000 to 0 

Got tho real data from block 
VP3I , channol VP32 and store it 
in VP2032 

Sot tho real data from block 
BP31 • channol VP33 and storo it 
in VP2033 

Got tho real data from block 
VP3I, channol VP3k and storo It 
In VP203A 

If VP2033<VP2030, skip 3 stops 

If VP2033<VP2032, skip 2 stops 

If VP2033<VP203k, skip I stop 

Sot tho complex data from block 
VP35, channol VP33 and storo it 
in VP 3000 

Put the complex data stored In 
VP3000 into block VP36, channol 
VP33 

VP32 • VP32 ♦ I 
VP33 • VP33 ♦ 1 
VP3<* - VP3k * I 

Loop through label 3k VP39 times 

VP35 - VP35 ♦ 2 

VP36 - VP36 + I 

VP3I • VP3I ♦ VP30 

Loop through label 33 three times 

Subreturn 

End 
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KEYBOARD FILE 4 - DIRECT POGO SOFTWARE 


1 L 

40 




U> 

X 

A 

18 




9 X> 

20 




13 Y _ 

3 

8 



19 Y _ 

4 

14 



25 Y 

8 

3 



31 Y _ 

10 

4 



37 Y _ 

12 

4 



43 Y _ 

14 

2 



49 Y _ 

16 

3 



55 Y _ 

18 

2 



61 L 

41 




65 CL 





68 CL 

1 




72 CL 

2 




76 CL 

3 




80 CL 

4 




84 CL 

5 




88 CL 

6 




92 Y 

41 

0 

0 

0 

99 Y A+ 

5 

-40 

1999 


106 Y A+ 

7 

-30 

1999 


113 L 

42 




117 Y 

45 

2 

0 

-3D 

124 Y 

45 

1 

0 

-4D 

131 * 

2 

-7D 



136 : 

2 

100 



141 * 

1 

-5D 



146 : 

1 

100 



151 SP 

1 

7 

2 



Label 40 

Load block 18 Into block 0 

Store block 0 Into block 20 

Set VP3 to 8 

Set VP4 to 14 

Set VP8 to 3 

Set VP10 to 4 

Set VP12 to 4 

Set VP14 to 2 

Set VP16 to 3 

Set VP18 to 2 

Label 41 


Clear blocks 0 through 6 


Call User Program 41, where 
nl»n2*n3«0 (Initialize off- 
line BSFA) 

VP5 • 1999 + VP(VP4) 

VP? « 1999 + VP(VP3) 

Label 42 

Call User Program 45; nl»2, n2“0, 
n 3«VP(VP3) (off-line BSFA) 

Call User Program 45; nl«l, 
n2«0, n3«VP(VP4) 

Multiply block 2 by VP(VP7) 

Divide block 2 by 100 

Multiply block 1 by VP(VP5) 

Divide block 1 by 100 

Calculate averaged 2-channel - 
double precision power spectrum 


50 


KEYBOARD FILE 4 - DIRECT 


157 # 

42 

2D 

163 CH 

1 

2 

169 X< 

1 


173 X> 

3D 


177 Y A+ 

3 


183 X< 

2 


187 X> 

3D 


191 : 

5 

2004D 

196 * 

5 

20 

201 : 

6 

2004D 

206 * 

6 

20 

211 X< 

5 


215 X> 

4D 


219 Y A+ 

4 


225 X< 

6 


229 X> 

4D 


233 Y A+ 

3 


239 Y A+ 

4 


245 # 

41 

3 

251 L 

43 


255 Y _ 

31 

9 

261 Y _ 

35 

8 

267 Y _ 

36 

5 

273 J 

32 

3 

278 L 

44 


282 Y _ 

3 

2 

288 Y _ 

4 

19 

294 Y _ 

5 

12 

300 Y _ 

6 

10 

306 Y _ 

7 

8 

312 Y _ 

8 

8 

318 Y 

9 

8 


B1I 

POGO SOFTWARE (Continued) 

Loop through label 42 VP2 times 

Calculate H(f) and y 2 (f) from 
the averaged power spectrum data 

Load block 1 Into block 0 [H(f)] 

Store block 0 Into block VP3 

VP3 * VP3 + 1 

Load block 2 Into block 0 [y 2 (f)3 

: ire block 0 in block VP3 

Divide block 5 by VP2004 

Multiply block 5 by 20 

Divide block 6 by VP2004 

Multiply block 6 by 20 

Load block 5 into block 0 

Store block 0 in block VP4 

VP4 = VP4 + 1 

Load block 6 into block 0 

Store block 0 into block VP4 

VP3 * VP3 + 1 

VP4 * VP4 + 1 

Loop through label 41 three times 

Label 43 

Set VP31 to 9 

Set VP35 to 8 

Set VP36 to 5 

Jump to label 33 in Keyboard File 3 

Label 44 

Set VPS to 2 

Set VP4 to 19 

Set VP5 to 12 

Set VP6 to 10 

Set VP7 to 8 

Set VP8 to 8 

Set VP9 to 8 


51 


BI2 


KEYBOARD FILE 4 - DIRECT POGO SOFTWARE (Concluded) 


324 Y __ 

10 

8 

Set VP10 to 3 

330 Y _ 

11 

15 

Set VP11 to 15 

336 Y _ 

12 

58 

Set VP12 to 58 

342 Y _ 

13 

150 

Set "PI 3 to 150 

348 J 

50 

5 

Jump to Table 50 in Keyboard 
rile 5 (Plot Routine) 

353 . 



End 


52 


KEYBOARD FILE 5 - POGO PLOT ROUTINE 


B1 3 


1 L 50 
5 Y 5838 3D 


10 

Y 


5821 

6 



15 

Y 


5865 

3 



20 

L 


51 




24 

Y 


26 

4D 



30 

Y 

IF 

4 

5D 

5 

0 

38 

Y 

IF 

4 

6D 

4 

0 

46 

Y 

IF 

4 

7D 

3 

0 

54 

Y 

IF 

4 

8D 

2 

0 

62 

Y 

IF 

4 

9D 

1 

0 

70 

Y 

IF 

4 

10D 

1 

2 

78 

J 


no 




82 

Y 


5809 

0 

0 

0 


89 

J 

100 



93 

J 

120 



97 

L 

52 



101 

Y A- 

4 



106 

# 

51 

1 ID 

0 

112 

L 

53 



116 

J 

130 



120 

Y 

5814 



124 

MS 

34 

1 2D 


129 

MS 

14 



133 

Y R 

4 




Label 50 

Call User Program 5838; nl=VP3 
(Read text buffer VP3 Into core) 

Call User Program 5821; nl»6 
(Set up the terminal as the 
plot device) 

Call User Program 5865; nl =3 
(Set log plots to plot 3 decades, 
vertical scale) 

Label 51 

Set VP26 to the value of VP4 

If VP4 * VP5, skip five steps 

If VP4 = VP6, skip four steps 

If VP4 * VP7, skip three steps 

If VP4 = VP8, skip two steps 

If VP4 = VP9, skip one step 

If VP4 > VP10, skip one step 

Jump to label 110 

Call User Program 5809; nl=0, 
n2=0, n3*0 (Set up to plot real/ 
magnitude, linear horizontal 
scale, scale switch at 12 
o'clock) 

Jump to label 100 
Jump to label 120 
Label 52 
VP4 = VP4 - 1 

Loop through label 51 VP11 times 
Label 53 

Jump to label 130 

Call User Program 5814 (Blank 
the terminal screen) 

Set the text record pointer to 
VPl 2 

Read the next message (plot 
menu) 

Input VP4 (data block number to 
be plotted 
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KEYBOARD FILE 5 - POGO PLOT ROUTINE (Continued) 


B14 


138 Y 


26 

4D 



Set VP26 to the value of VP4 

144 Y 

IF 

4 

10D 

2 

1 

If VP4 > VP10, skip 2 steps 

152 Y 

IF 

30 

0 

1 

0 

If VP30 = 0, skip 1 step 

160 Y 

A+ 

26 

26D 

30 


VP26 = VP26 + 30 

167 MS 

14 




Read the next message (Set the 
scale and press ‘■CONTINUE* 1 ) 

171 D 


4D 




Display block VP4 

175 Y 


£809 




Call User Program 5809; nl*n2* 
n3=default (Set up the plot 
according to the display 
switches) 

179 J 


100 




Jump to label 100 

183 J 


120 




Jump to label 120 

187 D 


4D 




Display block VP4 

191 J 


53 




Jump to label 53 

195 L 


100 




Label 100 

199 Y 


5814 




Call User Program 5814 

203 Y 


5829 




Call User Program 5829 (Plot 
the full data block) 

207 Y 


5815 

40 



Call User Program 5815; nl*VP4 
(Plot the data from block VP4) 

212 Y 


5816 




Call User Program 5816 (Plot 
the defaulted axis) 

216 Y 


5817 

26D 



Call User Program 5817; nl=VP26 
(Label the vertical axis with 
text buffer message VP4 and the 
horizontal axis with the 
default label) 

221 Y 


5808 

975 600 



Call User Program 5808; nl=970, 
n2=600 (Position the cursor to 
y*97.0" of the vertical plot 
window and x=60.0* of the 
horizontal plot window) 

227 Y 


5819 

1 



Write text buffer message #1 

232 < 






Subreturn 

235 L 


no 




Label 110 

239 Y 

IF 

4 

10D 

2 

1 

If VP4 >_ VP10, skip 2 steps 

247 Y 

IF 

30 

0 

1 

0 

If VP30 = 0, skip 1 step 


5 1 * 




KEYBOARD 

FILE ! 

> - 

POGO PLOT ROUTINE (Concluded) B15 

255 Y A+ 

26 

26D 

30 


VP26 = VP26 + 30 

262 TP 

4D 




Convert block VP4 to polar 
coordinates 

266 Y 

5809 

1 

0 

0 

Call User Program 5809; nl*l, 
n2=0, n3=0 (Set up to plot 
imaginary /phase, linear hori- 
zontal, scale switch to 12 
o'clock) 

273 J 

100 




Jump to label 100 

277 Y 

5808 

970 

1 3D 


Call User Program 5808; nl=970, 
n2=VPl 3 

283 Y 

5819 

260 



Call User Program 5819; nl=VP26 

288 J 

120 




Jump to label 120 

292 Y 

5809 

0 

0 

0 

Call User Program 5809; nl=0, 
n2=0, n3=0 

299 J 

130 




Jump to label 130 

303 J 

100 




Jump to label 100 

307 J 

120 




Jump to label 120 

311 TL 

40 




Convert block VP4 to log 
polar coordinates 

315 Y 

5809 

0 

0 

0 

Call User Program 5809; nl=0, 
n2=0, n3=0 

322 J 

130 




Jump to label 130 

326 J 

100 




Jump to label 100 

330 J 

120 




Jump to label 120 

334 J 

52 




Jump to label 52 

338 L 

120 




Label 120 

342 J 

130 




Jump to label 130 

346 Y 

5820 




Call User Program 5820 (Hard 
copy the terminal) 

350 < 





Subreturn 

353 L 

130 




Label 130 

357 Y _ 

19 

1 



Set VP1 9 to 1 

363 L 

131 




Label 131 

367 Y A+ 

19 




VP19 = YP19 + 1 

372 « 

131 

50 

0 


Loop through label 131 fifty times 

378 < 





Subreturn 

381 . 





End 


55 



Line * 


0 

1 

2 

3 

4 

5 

6 
/ 

8 

9 


10 

11 


12 

13 


14 

15 


20 

21 


APPENDIX C 

P0G0 SOFTWARE TEXT RECORDS (FILE 4a) 


Cl 


Contents 
Message #1 


This is a 5451C program, written for POGO data analysis. The program 
reads 4 analog channels simultaneously and writes them directly to the 
disc via the ADC Throughput. The time domain data is then called from 
the disc and H(f)s calculated re. Pulser. H(f)s across devices are 
calculated by ratios. 28 JANUARY 80 

Setup the ADC for the desired frequency resolution (Delta f), and 
input signal-to-noise ratio. Then, enter Blocksize. (INTEGER LE 1024) 

/* 


Message *2 

Enter a calibration factor for each channel. (4 - FLOATING POINT) 
/* 


Message *3 

Enter the number of data acquisition loops. (INTEGER) 

/* 


Message *4 

Enter Delta f as specified by the ADC. (FLOATING POINT) 

/* 

i 


Message *7 

Start the data tape. Press "CONTINUE." 

/* 


57 


C2 

Line * 


Contents 


Message ?8 

RATIO P060 DATA PLOT MENU 

24 The program output is located in core memory as follows: 

25 

26 BLOCK* DATA 

27 


28 

3 

Weighted H(f) HPOP/LPOP 

29 

4 

Weighted H(f) MCC/HPOP 

30 

5 

Weighted H(f) MCC/LPOP 

31 

6 

H(f) 

HPOP/LPOP 

32 

7 

Coh(f) 

ii 

33 

8 

H(f) ' 

MCC/HPOP 

34 

9 

Coh(f) 

it 

35 

10 

H(f) 

MCC/LPOP 

36 

11 

Coh(f) 

ti 

37 

12 

H(f) 

LPOP/Pulser 

38 

13 

Coh(f) 

H 

39 

14 

H(f ) 

HPOP/Pul ser 

40 

15 

Coh(f) 

it 

41 

16 

H(f) 

MCC/Pul ser 

42 

17 

Coh(f) 

M 

43 

18 

PSD 

Pulser 

44 

19 

PSD 

LPOP 

45 

20 

PSD 

HPOP 

46 

21 

PSD 

MCC 

47 

22 

Re XPSD LPOP/Pulser 

48 

23 

Im XPSD 

ii 

49 

24 

Re XPSD HPOP/Pulser 

50 

25 

Im XPSD 

it 

51 

26 

Re XPSD MCC/Pul ser 

52 

27 

Im XPSD 

ii 

53 




54 

Please 

enter the 1 

31 ock* you wish 


55 /* 

Message =9 

56 Set the scale and press "CONTINUE." 

57 /* 


58 


C3 


Line # Contents 


Message #10 

OIRECT POGO DATA PLOT MENU 


58 

The program output is located in core as fol 

59 



60 

BLOCK# 

DATA 

61 



62 

5 

Weighted H(f) HPOP/LPOP 

63 

6 

Weighted H(f) MCC/HPOP 

64 

7 

Weighted H(f) MCC/LPOP 

65 

8 

H(f ) HPOP/LPOP 

66 

9 

Coh(f) 

67 

10 

H(f) MCC/HPOP 

68 

11 

Coh(f) 

69 

12 

H(f) MCC/LPOP 

70 

13 

Coh(f) 

71 

14 

Re XPSD HPOP/LPOP 

72 

15 

Im XPSD 

73 

16 

Re XPSD MCC/HPOP 

74 

17 

Im XPSD 

75 

18 

Re XPSD MCC/LPOP 

76 

19 

Im XPSD 

77 



78 

Please enter the block # you wish to plot. 

79 

/* 




Message #11 

80 

Set the scale and press "CONTINUE". 

81 

/* 




Message *12 

82 

Do you want: Coherence blanking (enter 2) 

83 


Pulser blanking (enter 0)? 

84 

/* 



Message *13 

85 Enter the minimum blanking value (floating point). 

86 /* 


59 



APPENDIX D 


POGO SOFTWARE TEXT BUFFERS (FILE 4b) 


Text Buffer #1 
RATIO POGO PLOT LABELS 

01 

(TEST ID) 

/* 

13 

HPOP/LPOP Ratio H(f) 

Pulser Blanked 

/* 

04 

MCC/HPOP Ratio H(f) 

Pulser Blanked 

/* 

05 

MCC/LPOP Ratio H(f) 

Pulser Blanked 

/* 

06 

HPOP/LPOP Ratio H(f) - Polar 

/* 

07 

HPOP/LPOP Ratio Coh(f) 

/* 

08 

MCC/HPOP Ratio H(f) - Polar 

/* 

09 

MCC/HPOP Ratio Coh(f) 

/* 

10 

MCC/LPOP Ratio H(f) - Polar 

/* 

11 

MCC/LPOP Ratio Coh(f) 

/* 


61 


D2 

Text Buffer #1 (Continued) 

12 

LPOP/Pulser H(f) - Polar 

/* 

13 

LPOP/Pulser Coh(f) 

r 

14 

HPOP/Pulser H(f) - Polar 

/* 

15 

HPOP/Pulser Coh(f) 

/* 

16 

MCC/Pulser H(f) - Polar 

r 

17 

MCC/Pulser Coh(f) 

r 

18 

Pulser PSD 

/* 

19 

LPOP PSD 

r 

20 

HPGP PSD 

/* 

21 

MCC PSD 

/* 

22 

LPOP/Pulser Real XPSD 

/* 

23 

LPOP/Pulser Imag XPSD 

r 

24 

HPOP/Pulser Real XPSD 

/* 


62 


03 


Text Buffer #1 (Concluded) 

25 

HPOP/Pulser Imag XPSD 

/* 

26 

HCC/Pulser Real XPSO 

/* 

27 

HCC/Pulser Imag XPSO 

/* 

33 

HPOP/LPOP Ratio H(f) 

Coherence Blanked 

/* 

34 

MCC/HPOP Ratio H(f) 

Coherence Blanked 

/* 

35 

MCC/LPOP Ratio H(f) 

Coherence Blanked 

/* 


Text Buffer #2 
DIRECT POGO PLOT LABELS 

01 

(TEST ID) 

/* 

05 

HPOP/LPOP Direct H(f) 

Pulser Blanked 

/* 

06 

MCC/HPOP Direct H(f) 

Pulser Blanked 

r 

07 

MCC/LPOP Direct H(f) 

Pulser Blanked 

/* 
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Text Buffer #2 (Continued) 


Ok 


08 

HPOP/LPOP Direct H(f) - Polar 

/* 

09 

HPOP/LPOP Direct Coh(f) 

r 

10 

MCC/HPOP Direct H(f) - Polar 

r 

n 

MCC/HPOP Direct Coh(f) 

/* 

12 

MCC/LPOP Direct H(f) - Polar 

/* 

13 

MCC/LPOP Direct Coh(f) 

/* 

14 

HPOP/LPCP Direct Real XPSD 

r 

15 

HPOP/LPOP Direct Imag XPSD 

/* 

16 

MCC/HPOP Direct Real XPSD 

r 

17 

MCC/HPOP Direct Imag XPSD 

r 

18 

MCC/LPOP Direct Real XPSD 

r 

19 

MCC/LPOP Direct Imag XPSD 

r 

35 

HPOP/LPOP Direct H(f) 
Coherence Blanked 


64 


1 


Text Buffer #2 (Concluded) 

36 

MCC/HPOP Direct H(f) 

Coherence Blanked 

/* 

37 

MCC/LPOP Direct H(f) 

Coherence Blanked 


i 
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